Quantum light depolarization: the phase-space perspective 
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Quantum light depolarization is handled through a master equation obtained by coupling dispersively the field 
to a randomly distributed atomic reservoir. This master equation is solved by transforming it into a quasiproba- 
bility distribution in phase space and the quasiclassical limit is investigated. 
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I. INTRODUCTION 

Polarization of light is a key concept that has deserved a lot 
of attention over the years. Apart from its fundamental sig- 
nificance, it is also of interest in several active technological 
fields. In many of these applications it is crucial to determine 
the decrease of the degree of polarization experienced by the 
light when traversing an optical system: we refer to this as 
depolarization [j]]. 

In classical optics, this depolarization is ascribed, broadly 
speaking, either to birefringence (as it usually happens, e.g., in 



optical fibers Hi, HJ|, 
distributed particles 



]) or to scattering by randomly 
1 111 . In both cases the net result 
is an effective anisotropy that leads to a decorrelation of the 
phases of the electric field vector. 

In quantum optics, a sensible approach to deal with this 
decorrelation is through the notion of decoherence, by which 
we loosely understand the appearance of irreversible and un- 
controllable quantum correlations when a system interacts 
with its environment fl2tl . Usually, decoherence is accom- 
panied by dissipation, i.e., a net exchange of energy with the 
environment. However, we are interested in the case of pure 
decoherence (also known as dephasing), for which the pro- 
cess of energy dissipation is negligible. Models in which the 
populations do not change, while the coherences are strongly 
decaying, are at hand lfl3i [Til [l5tl . Since a good knowledge 
of dephasing is of utmost importance (prominently for quan- 
tum information processing, where operations completely rely 
on the presence of coherence), these models have been suc- 
cessfully applied to dephasing in systems such as quantum 
dots CI O [H [Hi, Josephson junctions HJ HE [22, [H, 
or general quantum registers l24ll . to cite only a few relevant 
examples. 

In the modern parlance of quantum information it is usual to 
call depolarizing channel a decoherence induced by an unbi- 
ased noise generating bit-flip and phase-flip errors [25]. While 
this terminology fits well with abstract qubits, whenever light 
is concerned there is an extra essential ingredient to be taken 
into account: quantum polarization has a quite natural su(2) 
invariance that leads to a structure of invariant subspaces 12611 . 
On physical grounds, we argue that this structure must be pre- 
served in the evolution, which makes previous approaches to 
fail in this case. 

The purpose of this paper is twofold. First, we wish to pro- 



vide a simple model that goes around this drawback and gives 
a picture of the mechanisms involved in depolarization. The 
main idea is to couple the field dispersively to a randomly dis- 
tributed atomic reservoir: the resulting master equation has a 
quite appealing structure that complies with our requirements 
of su(2) invariance. Our second goal is to solve this equation 
in phase space, presenting then a representation of the depo- 
larizing dynamics that makes an easy contact with the classi- 
cal one on the Poincare sphere. 



II. SU(2) POLARIZATION STRUCTURE 

We assume a monochromatic plane wave propagating in the 
z direction, whose electric field lies in the xy plane. Under 
these conditions, we have a two-mode field that can be de- 
scribed by two complex amplitude operators. They are de- 
noted by a+ and <z_, when using the basis of circular (right 
and left) polarizations, which we shall employ in this paper. 
The commutation relations of these operators are standard: 



[a x ,a\,] = 5\x>i, A, A' e {+,-}. 



(2.1) 



The Stokes operators are then defined as the quantum coun- 
terparts of the classical variables, namely JZ7il28[|29l[30l[3Tll 



So = a} + a + + ala_ 
§2 = i(a+a_ — a\_d- 



S\ = a\ a_ + ciLa + , 



(2.2) 



S, 



it 



ci- 



and their mean values are precisely the Stokes parameters 
((S ), (S», where S = (S\ , S 2 , S 3 ). Using the relation dXTI 
one immediately gets that the Stokes operators satisfy the 
commutation relations 



[S,S ] = 0, [S 1 ,S 2 ] = 2iS 3 



(2.3) 



and cyclic permutations. Since So is just the operator repre- 
senting the total number of photon, the first equation in ( 12.31 ) 
means that S is measurable in photon-counting experiments 
and we can treat each energy manifold separately. To bring 
out this point more clearly, it is advantageous to relabel the 
standard two-mode Fock basis in the form 



\N, k) = \k)+® \N-k)- 



0,1, 



,N. (2.4) 



2 



For each fixed value of the number of photons N, these states 
span an invariant subspace of dimension N + 1 and the oper- 
ators S act therein according to 

S+\N,k) = 2y/(k + l)(N -k)\N,k + l), 
S-\N,k) = 2yfk{N - k+ T)\N,k- 1) , (2.5) 
S 3 \N,k) = 2{k-N/2)\N,k), 

where S± = Si ± iS2 are raising and lowering operators. 
These invariant subspaces will play a key role in the following. 

The quantities S are the generators of the su(2) algebra. The 
noncommutability of these operators precludes the simultane- 
ous exact measurement of their physical quantities. Among 
other consequences, this implies that no field state (leaving 
aside the two-mode vacuum) can have definite nonfluctuating 
values of all the Stokes operators simultaneously. This is ex- 
pressed by the uncertainty relation 

(AS) 2 = (AS*i) 2 + (AS 2 ) 2 + (A5 3 ) 2 > 2(S ) . (2.6) 

Contrary to what happens in classical optics, the electric vec- 
tor of a monochromatic quantum field never describes a defi- 
nite ellipse. 

Finally, we recall that standard definition of the degree of 
polarization reads as 

m |(S)| (Si) 2 + (S 2 ) 2 + (S 3 ) 2 

P = = * . (2.7) 

(So) {So) 

Note that P depends exclusively on the first moments of the 
Stokes operators. Higher order moments can be crucial for a 
full understanding of some phenomena ll32il33i[34ll . but ( 12.7b 
is more than enough for our purposes here. 

III. DEPOLARIZING DYNAMICS IN A NONRESONANT 
ATOMIC MEDIUM 

Pure dephasing dynamics occurs when the system Hamil- 
tonian _ff sys commutes (at least, approximately) with the 
system-reservoir interaction or for initial states that evolve 
very slowly under the dynamics governed by -ff sys on the time 
scale of decoherence processes 13511 . Various theoretical sce- 
narios have been proposed to deal with this: apart from minor 
details, all of them can be modelled by a scattering process in 
which a reservoir quantum can be absorbed or emitted, but the 
number of excitations in the system is unchanged [13]. These 
models indeed preserve the invariant subspaces, but produce 
no thermalization: they merely maintain the occupation prob- 
abilities, while erasing all coherences. Therefore, they fail to 
describe light depolarization because depolarization not only 
preserves the invariant subspaces, but the steady state in each 
one of them must be a diagonal state (36J[. We believe that 
these conditions are essential to ensure a correct description 
of any depolarizing process. 

To solve these difficulties we take another route: we assume 
that the field mode of frequency to propagates through a ma- 
terial medium made of a collection of two-level atoms with 



transition frequencies uj a . In consequence, our basic system 
is represented by (in units h = 1 that we shall use throughout 
all this paper) 

H sys = Hf + H a + Hf a , (3.1) 

where 

X=± 
a 

= (- 9aA °« 4 + 9* a a °t s a) • 

a A=± 

Here, as is normal in practice, we have described each indi- 
vidual atom in terms of the standard Pauli operators. The 
sum over a runs over all the atoms and the interaction H{ a 
is written in the dipolar and rotating-wave approximations. 
In addition, the atoms are randomly distributed so the cou- 
pling constants g a \ have random phases: we thus write them 
as g a \ = \Qa\e iXVa ^ 2 , for A = ±. This random-phase ap- 
proximation 03711 is used in many areas of physics to esti- 
mate response functions and it works properly in the long- 
wavelength limit, a hypothesis implicit in the form of the 
Hamiltonian (|3.2| i. 

It is a well-known fact that atoms decay irreversibly. This 
is usually assigned to their interaction with the continuum of 
modes of an additional thermal electromagnetic environment. 
When we take this point into account, the density matrix for 
the system d3.2| > evolves according to 

Q S yM = -i[H sys , g sys ] 

+ ~2~{( na + &ys + n a C[a+] g sys } , 

a 

(3.3) 

where C[C] are Lynblad superoperators J38ll 

C[C] g = 2Cq& -{C' t C,g}, (3.4) 

and 7 a is the decay constant of the ath atom due to its coupling 
to the thermal environment with h a excitations. Note that, 
by introducing a different bosonic reservoir for each atom, 
we avoid the occurrence of any collective effect. Whereas 
a collective dephasing can introduce remarkable dynamical 
changes JH,|43,|4ll|43,|4l|4l], we expect it plays no relevant 
role in understanding the quasiclassical limit of depolarizing 
dynamics we are interested here ll45ll . 

The properties of a random medium are well reproduced 
when n a ^> 1 j46l l47ll . In this high-temperature limit the 
effect of spontaneous emission can be disregarded in compar- 
ison with stimulated processes. Emission into the reservoir 
and absorption from the reservoir therefore become identi- 
cal; i. e., they balance each other in the stationary state: the 
emission and absorption processes depend solely on the initial 
population of the atomic state. Consequently, the steady-state 
reduced density operator is approximately given by a mixture 



3 



of equally populated atomic states and the density matrix of 
the ath atom becomes thus diagonal (g a — it). 

If A Q = u> a — u> is the detuning, we consider the far off- 
resonant regime 



« |AJ 



(3.5) 



In this limit we can adiabatically eliminate the atomic vari- 
ables and obtain a master equation that, after averaging over 
the random phases, reads as (see Ref. l48tl for technical details 
on the derivation) 

'Q = -iuj[S ,Q]+2jjC[§ }g+-fjC[S + ]g+jC[^}g, (3.6) 

where g(t) — Tr at [p sys (i)] is the reduced density operator for 
the field mode and the decoherence rate 7 is 



7 = E 



la A 2 n a 



(3.7) 



The master equation ( 13.6b preserves the N -photon subspaces 
and the steady state in each one of them is a completely ran- 
dom state 



g(t — > 00) 



1 



N + l 



11 , 



(3.8) 



which is an important benefit of our approach. The depolar- 
ization rates ( 13.71 i are very small when compared with other 
typical system parameters (observe the dependence on A~ 2 ), 
in agreement with experimental observations. 

The terms £[S±] describe depolarization in each invariant 
subspace, meanwhile the action of C[Sq] therein is trivial. 
Nevertheless, C[Sq] is responsible for the relative phase decay 
between blocks of the density matrix corresponding to differ- 
ent excitation numbers. 

As an illustrative example, we focus on one-photon states. 
We are then within a two-dimensional subspace and the den- 
sity matrix can be expressed as 



Q 



If 1 + z x — iy 
2\x + iy 1 — z 



1 



(11 



(3.9) 



where r = 2 Tr (gcr), cr being the Pauli matrices. The vector 
r satisfies |r| < 1, with the equality for pure states (that are 
all them SU(2) coherent states ll49ll ). 

In this subspace the master equation ( 13.6b can be recast in 
terms of r(t) and one obtains the solution as 

x{t) = a;(0)e- 7 *, 

y(t) = y(p)e-T, (3.10) 
z(t) = z(0)e- 27t . 

If we apply the definition ( 12.7b to this case, we get 



-|r(t)|e-T*. (3.11) 



This degree tends then to zero with a typical time scale 7 1 , 
which, as commented before, is exceedingly large. 



IV. THE QUASICLASSICAL REGIME 

The first term in the master equation (13.6b produces only 
free evolution and can be dropped in what follows. In conse- 
quence, we consider the depolarizing equation 



g = TC[S Q ] g + 7 £[S+] g + 7 £[5_] g , 



(4.1) 



where, for the sake of generality, we have considered different 
decoherence rates for So and S± . To work out the solutions of 
(14. lb we resort to phase-space methods, since they provide the 
most suitable approach to study the quasiclassical limit. Ac- 
cording to the discussion in Sect. [TU one could think in using 
quasiprobability distributions on the sphere f50l l5ll HH l53ll . 
However, although naturally related to the SU(2) group, they 
are not so widely employed by the quantum-optics commu- 
nity, so we prefer to take a more standard approach. There- 
fore, we shall use the s-parametrized quasiprobability distri- 
butions introduced by Cahill and Glauber 15411 . which for our 
two-mode fields read as 



W«(o+,a_) = i J d 2 /3 + d 2 /3_x (s) (/? + ,/3-) 

^2 - a *\P> 



x exp 



(4.2) 



where ) is the s-ordered characteristic function 

X (s) (/?+ , P- ) - Tr [g (4.3) 
and £>( s ) (J3+, is the two-mode displacement operator 
D^(P+,P-) = J] e s ^ 2 / 2 exp(j3 x a{ - fcax) . (4.4) 



In the three special cases s = 0, +1, and —1, one easily 
recognizes the Husimi Q, the Wigner W, and the Glauber- 
Sudarshan P functions, respectively. In terms of W^ s \ 
Eq. (14.1b can be transformed into a differential equation us- 
ing the following rules of mapping J55ll 



a\g 



a[g 



ga x 



ga{ 



1 




s 


d 




2 






1 


+ 


s 


d 




2 




da\ 


1 


+ 


s 


d 




2 




da* x 


1 




s 


d 




2 




da x 



w {s) 
w {s) 



(4.5) 



For definiteness, we choose to work with the P function, 
which will facilitate the following calculations. Applying the 
rules J4.5b for each mode, we get that Eq. J4.ll ) can be equiv- 
alently written as 



4 



P(a+,ct-) = — r I a+ d a+ — a* + d a * — a.- d a _ + a*_ d a *_ ) P(a+,a_) 



2" 



\a-\ 2 d a+ d a * + + \a + \ 2 d a _d a *_ - a+a- d a+ d a _ - a* + a*_ d a * + d a * 



i (a+ d a+ + a* + d a - + + a_ d a _ + a* d a * 



P(a + ,a-) 



(4.6) 



r 



At first glance, this equation looks very intricate to be of any 
practical value. However, let us introduce the following dif- 
ferential operators 



S 


= a- 


d a _ 


— a_ 


d a *_ 


+ a + d a+ 


s± 


= a T 


d a± 


* 

-a ± 






s 3 


= a_ 


d a _ 


-a*_ 




- a + d a+ 



If we concentrate in the evolution of the components of the 
Bloch vector (S), we have that 

(a + ,a^\S 3 \a + ,a^) = r 2 cos ^ r 2 Dl ((f), 9, ip) , 

(4.12) 

(4.7) (a+,a_|£±|a + ,a_) = r 2 e T ^ sin<9 = \/2r 2 D^ty, 6, ijj) , 



a*+ d a * 



which are nothing but a differential realization of the Stokes 
operators (i.e., of the u(2) algebra). It turns out that Eq. ( 14.6b 
reduces to the very simple factorized form 



P(a+,a-) 



T5 2 + ~f(S+S- + S-S+) P(a+, a 



(4.8) 

But S 2 = <S| + (<S+<S_ +S-S + )/2 is the Casimir operator for 
su(2), which reduces to 5(5 + 1) 1 in each invariant subspace 
(5 taking only integer or half-odd integer values). Then, this 
master equation has the formal solution 

P(a+,a-\t) = cxp{-4[r5 2 + 2 7 (5 2 -4 2 M 

x P(a + ,a-\t = 0). (4.9) 

Because of the properties of the P representation, the average 
value of any observable 0(t) can be recast as 

(6(t)) = J d 2 a+d 2 a-P(a +7 a-\t = 0) 

x exp{-4[r5^ + 2 7 (5 2 - S^] (a+a^ \6\a+aJ) , 

(4.10) 

where by O we mean the operator at t = 0. The problem for 
a practical application of this equation is that we need to give 
meaning to the action of the exponential. To this end, we first 
introduce the parametrization 



a+ = re- 4( ^ )/2 cos(e»/2) 



(4.11) 



where r 2 = |a+| 2 + |ct-| 2 is a radial variable related with 
the global intensity. The angles 8 and tf> can be interpreted as 
the polar and azimuthal angles, respectively, on the Poincare 
sphere, while -ip is a common phase. These coordinates reflect 
the fact that polarization needs only three independent quan- 
tities to be fully characterized: the amplitudes of each mode 
and the relative phase between them. 



where Df nm , are Wigner /^-functions IdoI . which constitute 
an orthogonal basis of complex-valued functions defined on 
the irreducible carrier subspaces of su(2), so that the action of 
the operators (14.7b on them is standard. Then a direct applica- 
tion of (14.101 1 yields 



(S 3 (*)> = exp(-8 7 t)(5 3 (0)) 



(S±(t)) = exp[-(27 + r)t](S±(0)>. 



(4.13) 



The procedure can be applied much in the same way for the 
field amplitudes; we merely quote the simplest results: 



(a±(t)) = exp[-(7 + r/4)t]<a ± (o)>, 



(a 2 ± (t)) = exp[-2( 7 + r)t](a 2 fc (0)) 



(4.14) 



The rest of the moments can be obtained in similar fashion. 

On the other hand, in terms of ( 14.1 It the operators So and 
£3 take the suggestive form l56ll 

So = -idip , S 3 = -id,/, , (4.15) 
and the exact solution of ( |4.8t can be expressed as 

00 s 

p^,e,m = E E ^C^-w) 

S m,m' — — S 

x exp{ [-87 5(5 + 1) + 8 7 m 2 - 4Tm' 2 ]t} , 

(4.16) 

where the coefficients cf n , are determined by the initial state: 

87T 2 /" 27r c K f 2rr 

4am' = 2S+1 J ^ J o d9 S[lld J q #Cm'WtM) 



X P((f),e,1p\t = 0) 



(4.17) 



We can proceed to integrate over the physically irrelevant 
global phase if>, obtaining in this way the P function over the 



5 



unit sphere. The result is 



p(<t>,6\t) = EV^TT ^ c s m0 Y Sm (o,4>) 

S m= — S 

x exp{-8 7 [S'(S'+l)-TO 2 ]t}, (4.18) 

where Yg m (9,(/)) are spherical harmonics. This is a closed 
formula that can be considered as our major result: it allows 
to trace the depolarizing dynamics on the Poincare sphere for 
arbitrary states. For example, for (quadrature) coherent states, 
we have 



coefficients c^ m , are determined now as in ( 14.17b . but with P 
replaced by Q] . 

This can be used to determine any function of Q. In fact, 
it has been argued ll57ll that, for a given field, the distance be- 
tween its Q function and the Q function for unpolarized light 
can be taken as a proper degree of polarization in phase space: 
this distance is, except from constant factors, proportional to 
Q 2 . Therefore, we can take as an unnormalized depolarization 



/•27T i>TT />27T 

Kt)= dip d0 sin. d<f>Q 2 (<j),$ 7 iP\t), (4.22) 
Jo Jo Jo 



P(a + ,a-\t — 0) = 5 2 (a + - r )S 2 (a- - r ) , (4.19) and then we obtain 



where we have assumed that both modes have the same real 
amplitude ro. We then get 

P(<P, 9, ip\t = 0) = 8V2 ~ r °) 5{i>) S(</>) 6(9 - tt/2) , 

(4.20) 



oo . S 



„S |2 



m,m' = — S 



x exp[-8(2 7 S'(5 + 1) - 2 7 m 2 + lW 2 )](4.23) 

Again we can ascertain the corresponding time evolution for 
arbitrary states. 

V. CONCLUDING REMARKS 

In this paper we have reported a comprehensive, simple the- 
ory of quantum light depolarization. In our model the field 
couples dispersively to a randomly distributed atomic reser- 
voir, and the resulting master equation has unique properties 
that we have explored in detail. 

We have solved this master equation resorting to a sim- 
ple phase-space formalism for polarization on the Poincare 
sphere, based on s-ordered quasidistributions for the two ba- 
sic polarization modes. These results may have interesting 
consequences to implement experimental procedures for de- 
termining polarization properties. 



and so 

j3 



64V2vr 2 
25+1 



2^,(0,71/2,0) 



8(r- r ) 



(4.21) 



From here we can immediately compute the dynamics of this 
state. 

In view of the delta functions appearing in (14.201 , one can 
argue that the P representation results in a distribution with 
singularities. We can look instead at smoother quasidistribu- 
tions, such as, e.g., the Q function. Curiously enough, when 
we repeat the calculations, we end up with the fact that the 
evolution equation for the Q function is exactly the same as 
for the P function. In consequence, the time evolution of the 
Q function can be expressed as in Eq. (14.161 1 [obviously, the 
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